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Lensfree on-chip imaging and sensing platforms provide compact and cost-effective designs for various 
telemedicine and lab-on-a-chip applications. In this work, we demonstrate computational solutions for 
some of the challenges associated with (i) the use of broadband, partially-coherent illumination sources for 
on-chip holographic imaging, and (ii) multicolor detection for lensfree fluorescent on-chip microscopy. 
Specifically, we introduce spectral demultiplexing approaches that aim to digitally narrow the spectral 
content of broadband illumination sources (such as wide-band light emitting diodes or even sunlight) to 
improve spatial resolution in holographic on-chip microscopy. We also demonstrate the application of such 
spectral demultiplexing approaches for wide-field imaging of multicolor fluorescent objects on a chip. These 
computational approaches can be used to replace e.g., thin-film interference filters, gratings or other optical 
components used for spectral multiplexing/demultiplexing, which can form a desirable solution for 
cost-effective and compact wide-field microscopy and sensing needs on a chip. 

The optical microscope has been a very important tool for various fields, including engineering, medicine and 
sciences. Over the last decades, various breakthrough advances pushed the limits of optical microscopy in 
terms of resolution, contrast, signal-to-noise ratio (SNR), imaging speed, etc. opening up new applications 
for light microscopy. However, most of these advanced optical microscopy techniques are limited to advanced 
research and clinical settings due to their relatively high cost, complexity and weight. Computational microscopy 
tools 1 " 18 , on the other hand, offer more accessible optical imaging platforms by replacing costly and complex 
optical hardware with image reconstruction algorithms. Lensfree on-chip imaging and sensing platforms form an 
important example of this research theme. Besides reducing the cost and optical complexity, these lensfree 
microscopy techniques also aim to provide field-portability as well as compatibility with microfluidic devices 
by enabling wide field-of-view and high-throughput imaging without sacrificing resolution, making them highly 
desirable for various telemedicine 19 " 22 and lab-on-a-chip 5 ' 23 " 26 applications. 

Despite their design simplicity, on-chip imaging approaches also present some challenges that need to be 
addressed to improve the imaging performance, including the spatial resolution 1,27,28 . The first one of these 
challenges is a result of placing the objects very close to the detector array (i.e., on a chip) as illustrated in 
Fig. La. This short vertical distance (e.g., z 2 < 1 mm) is necessary to maintain a decent SNR for both lensfree 
holographic and fluorescent imaging on a chip and also allows all of the active area of the detector array to serve as 
the imaging field of view (FOV), typically achieving a FOV of for example 15-30 mm 2 and >l-20 cm 2 using 
CMOS and CCD image sensors, respectively. However, this small z 2 choice of on-chip microscopy tools also 
causes the pixel size of the detector to be a potential limiting factor toward achieving a high spatial resolution 
without the use of any lenses. As a result, on-chip holograms of objects suffer from under-sampling and aliasing 
artifacts as illustrated in Fig. 2. Under the assumption that the illumination is sufficiently narrow-band (e.g., ~3- 
5 nm), these spatial under- sampling artifacts can be well addressed by implementing pixel super resolution 29 " 31 
methods on lensfree geometries based on e.g., source 29,31,32 or object shifting 24,33 . Another challenge with on-chip 
imaging is due to the absence of lenses; i.e., one needs to reconstruct diffracted fields of the objects, which 
constitute (i) partially- coherent fields as in the case of holography and (ii) incoherent fields as in the case of 
fluorescent objects. Depending on the imaging geometry and coherence properties of the diffracted object fields, 
iterative phase retrieval 34 " 36 techniques or compressive decoding 37 algorithms 38,39 can be used to digitally recon- 
struct the images of the objects from their intensity measurements. These previous reconstruction 
approaches 1,22,26,40 assumed a single wavelength of operation, and the effective bandwidth of the illumination 
source or the fluorescent emitters was approximated to be narrow, i.e., quasi-monochromatic. This assumption, 
while valid in most cases, could introduce image artifacts and resolution limitations in on-chip imaging; for 
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Figure 1 | Lensfree holographic and fluorescent computational 
microscopy set-ups. (a) A schematic of the partially coherent, on-chip 
holographic microscopy setup, (b) A schematic and (c) a photograph of 
the lensfree, on-chip multicolor fluorescence microscopy platform. 

example in holographic on-chip microscopy a large bandwidth for 
the source (e.g., 20-30 nm) would result in a reduced temporal 
coherence length, which would also reduce the effective numerical 
aperture that can be synthesized after a pixel super-resolution tech- 
nique is applied. 

Compared to highly coherent laser sources, in our on-chip holo- 
graphic imaging geometry partially coherent light sources such as 



light emitting diodes (LEDs) 41 are preferred to mitigate speckle and 
multiple reflection interference noise terms. Therefore, despite their 
limited temporal coherence, LEDs are highly desired for field-port- 
able pixel-super resolution microscopy 24,31 ' 32,42 . In addition to signifi- 
cantly reduced speckle and multiple reflection interference noise, 
such broader band sources also improve the light efficiency and allow 
shorter image acquisition times as well as longer object-to-sensor 
distances. Therefore, computational demultiplexing methods are 
highly desirable toward digitally relaxing the temporal coherence 
limitation of broadband sources in on-chip holographic imaging. 
Another area where such computational spectral demultiplexing 
approaches can be particularly useful is lensfree on-chip fluorescent 
imaging, where the incoherent fields of various fluorescent emitters 
of different colors can be resolved despite their spatial overlap at the 
sensor plane. 

Toward this end, here we demonstrate a set of computational 
spectral demultiplexing approaches for broad-band holographic 
and multicolor fluorescent on-chip imaging. The first of these tech- 
niques will focus on holographic on-chip microscopy, where we 
report a polychromatic phase retrieval approach to demultiplex the 
superimposed spectral components of the object holograms recorded 
with a broadband LED illumination. Based on a similar idea, we also 
demonstrate a field-portable holographic microscope performing 
multispectral holographic on-chip imaging with sunlight 43 illumina- 
tion over a wide field-of-view of e.g., 15-30 mm 2 . In this second 
approach, the superimposed holograms are digitally recovered with 
the knowledge of the spectral response curves of both the pixels of the 
CMOS sensor chip and the illumination spectra through the use of a 
compressive decoding method, which provides a decent fit since this 
inverse problem is underdetermined and the multispectral spatial 
information of our objects is naturally redundant. Finally, we report 
an on-chip multicolor fluorescent microscopy method with a wide 
imaging field-of-view of e.g., —2.4 cm 2 . In this fluorescent imaging 
setup, the incoherent signatures from the objects are detected by a 
color (i.e., RGB) sensor-array that is placed at < 0.5 mm away from 
the sample plane. Diffraction correction and spectral demultiplexing 
of these multicolor lensfree fluorescent signatures are simultaneously 
performed through a sparse signal recovery algorithm with the 
knowledge of the system point spread function (PSF) at different 
wavelengths and the sensor's spectral response curve. These compu- 
tational demultiplexing approaches can in general be useful to 
replace various color filters, gratings or other optical components, 
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Figure 2 | Effects of pixel size and illumination bandwidth on lensfree image reconstruction quality, (a and b) Effects of under- sampling and (a and c) 
lack of temporal coherence for lensfree holographic on-chip microscopy, (c and d) The use of pixel super resolution to overcome spatial under- sampling. 
The cross-sections along the first diffraction orders together with the recovery results demonstrate the effects of spatial sampling rate and spectral 
bandwidth on image quality. 
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providing a desirable solution for cost-effective and compact on-chip 
holographic and fluorescent microscopy even in field-conditions. 

Results 

When the illumination source for lensfree holographic on-chip 
microscopy is not monochromatic, as in the case of broadband 
LEDs, the recorded lensfree hologram (I to tai) can be represented as 
an incoherent superposition of weighted monochromatic holograms 
(I x .)> such that: 

Itotai = Ci I Xl + c 2 h 2 + c 3 I h . . . + q I x . , (1) 

where I x = \R Xi +S x f = \R Xi \ 2 + \S Xi \ 2 + 2Re{R Xi Sl}> R Xi is the ref- 
erence wave and S Xi is the object wave for each wavelength. The 
overall contribution (q) of each wavelength within the broad illu- 
mination spectrum is determined by (/) the specific spectrum of the 
illumination source and (it) the spectral responsivity of the sensor 
chip. This weighted summation of these monochromatic holograms 
causes a spatial smearing effect in the measured broadband holo- 
gram, decreasing the fringe visibility and causing some of the high 
resolution information embedded within these fringes to be buried 
under the background noise as illustrated in Fig. 2. Another way to 
explain this effect is that the temporal coherence length shrinks as the 
source bandwidth (AX) broadens, i.e., 

X 2 

Coherence length, L c ~ cen ^ r (2) 



where n is the refractive index of the medium and ^ cen ter is the centre 
wavelength of the source. This reduced temporal coherence length of 
a broadband source limits the maximum spatial resolution and 
causes reconstruction artifacts, which become more pronounced 
especially for sub-micron spatial features of an object. To address 
this important challenge, here we use a polychromatic phase retrieval 
(PolyPR) approach to digitally relax the temporal coherence limita- 
tion and reduce the spatial artefacts due to monochromatic treat- 
ment of the holograms created with broadband illumination. 

To demonstrate how this polychromatic phase recovery approach 
can significantly improve the reconstruction quality under wide- 
band LED illumination, we designed a holographic on-chip imaging 
experiment (as in Fig. la) where the spatial resolution is limited by 
the temporal coherence of the source. For this end, the lensfree 
holograms of a 1.5 um period phase grating that is illuminated with 
a 15 nm bandwidth LED are recorded with a CMOS sensor-array 
(1.12 um pixel pitch) and are used to create a single, pixel super- 
resolved hologram 24 ' 31 ' 32 . This resulting hologram is then processed/ 
reconstructed using (i) the regular phase recovery approach 22 under 
the assumption that the source is quasi-monochromatic, and (ii) the 
PolyPR algorithm, which is detailed in Fig. 3 (top), Supp. Fig. 2 and 
the Methods section. The comparison of the reconstructed images in 
each case is demonstrated in Fig. 3 (bottom), where PolyPR approach 
presents significant improvements over the regular phase retrieval 
results, especially for the phase images. 

We should emphasize that the recovery of a spectral image 
cube from a single broadband measurement is in general an 
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Figure 3 | Polychromatic Phase Retrieval (PolyPR). PolyPR improves the spatial resolution by relaxing the limitations due to broadband illumination 
and lack of sufficient temporal coherence. (Top) A block diagram explaining the PolyPR method. (Bottom) Regular and polychromatic phase retrieval 
results for a 1.5 |am-period grating, which is illuminated with wide (15 nm) and narrow (3 nm) bandwidth (BW) LEDs. 
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underdetermined problem. As a result, despite being superior to the 
broadband illumination followed by regular PR case, PolyPR based 
results cannot match the performance of the narrowband illumina- 
tion under regular PR case. Also, PolyPR requires more iterations to 
converge compared to regular PR (e.g., —100 vs. 15), which may 
make some of the edge artefacts of object support-based phase 
retrieval technique more pronounced. This, however, can be miti- 
gated by using different constraints such as multi-height intensity 
measurements, which do not need any object support information 
for phase retrieval 42 . 

After these initial imaging results obtained with a broadband LED, 
we tested our technique by using the sun as our illumination source. 
Sunlight had been frequently used as a source in optics experiments 
until other light sources such as incandescent light bulbs were 
invented. Providing an easy-to-access and rather broad spectral con- 
tent, the sunlight still presents a rich source 11 that could potentially be 
useful in multispectral imaging to extract the wavelength- dependent 
features of microscopic objects of interest. In this context, to invest- 
igate the use of sunlight as an illumination source for our on-chip 
holography set-up, we created a portable holographic microscopy 
platform working with sunlight (see Fig. 4 and the Methods section 
for details). In this hand-held platform, sunlight based lensfree holo- 
grams are separated into multiple narrowband holograms by using a 
compressive sampling/sensing based algorithm, which is summar- 
ized in Fig. 5 (top) and is further detailed in the Methods section. The 
experimental performance of this platform is demonstrated by 
imaging a US Air Force Resolution Target. Our spectral demultiplex- 
ing approach clearly improved the quality of the recovered images 
(see Fig. 5, bottom), where we can reach a spatial resolution of 
—2.3 um over a field of view of —15 mm 2 using our field-portable 
holographic sunlight microscope (Fig. 4). 

In addition to their use for lensfree holographic on-chip imaging, 
as discussed in our Introduction, computational spectral demulti- 
plexing methods can also play a critical role for improving the 
performance of multicolor on-chip fluorescent microscopy. Per- 
forming fluorescence imaging without lenses by using an on-chip 
geometry provides a very large field- of- view, easily reaching e.g., 
> 10-20 cm 2 25 ' 26 . To further increase the multiplexing capabilities 
and the specificity of on-chip fluorescent imaging, performing 



simultaneous multicolor detection is highly desirable. For this pur- 
pose, we created a lensfree fluorescent on-chip microscopy plat- 
form, as shown in Figs. Lb and c, to image multicolor fluorescent 
emitters simultaneously. We also developed a decoding method 
(summarized in Fig. 6, with further details given in the Methods 
section) that simultaneously recovers the position and the color of 
the fluorescent emitters. The performance of this on-chip fluor- 
escent imaging approach is tested with 4 um green (excitation/ 
emission: 505/515 nm) fluorescent beads and 10 um red (excita- 
tion/emission: 580/605 nm) and green (excitation/emission: 505/ 
515 nm) fluorescent beads as shown in Fig. 7. The raw fluorescent 
signatures were modulated with the Bayer pattern and exhibit severe 
spatial overlap due to the lensfree nature of our on-chip imaging 
platform (see e.g., Fig. 7.al, bl, c5 and d5). These multicolor fluor- 
escent objects and their spatial distribution can be successfully 
recovered as illustrated in Fig. 7, which is cross-validated using 
fluorescent and bright field microscope images of the same objects. 
This simultaneous multicolor reconstruction feature further boosts 
up the throughput and the multiplexing capability of our wide-field 
on-chip imaging platform, which can be useful for various applica- 
tions including for example imaging and screening of rare cells 
using large volume micro -fluidic chips. 

Discussion 

Some of the limitations of our on-chip holographic imaging plat- 
forms due to a broadband illumination source could be mitigated or 
at least partially reduced by modifications of the optical setup such as 
having a shorter illumination wavelength (Center )> a closer object- 
to-sensor distance (z 2 )> or using spectral filters to narrow the illu- 
mination bandwidth. However, these physical modifications to the 
on-chip holographic imaging set-up also have their own short- 
comings. First, most image sensors exhibit reduced sensitivities at 
shorter wavelengths. Second, reducing the object-to-sensor distance 
causes severe under- sampling of the holographic signatures, and 
might also yield insufficient lateral shifts for pixel super resolution 
techniques 3 . Another potential issue for the use of smaller object-to- 
sensor distances is reduced thermal isolation from the sensor- array, 
which might create limitations for imaging of e.g., live biological 




Figure 4 | A field-portable, on-chip platform is designed to perform holographic imaging with sunlight, (left) A detailed schematic and (right) a photo 
of our portable, on-chip platform for holographic imaging using sunlight. 
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Figure 5 | Spectral demultiplexing of lensfree sunlight holograms using compressive decoding. (Top) A block diagram explaining our decoding method 
used to demultiplex the lensfree sunlight holograms into narrower spectral bands. (Bottom) Back-projection results are reported for lensfree holograms 
of a resolution test chart, before and after the spectral demultiplexing step. The inset images clearly show the improvement, (bottom, right inset) 
where a spatial resolution of —2.3 (am is demonstrated. 



objects. Therefore, digital narrowing/demultiplexing of the band- 
width of the source is a rather convenient, cost-effective and 
photon -efficient solution that can be used to partially overcome 
the resolution losses due to the limited temporal coherence of broad- 
band or polychromatic sources in lensfree on-chip holographic 
microscopy. To demonstrate this idea, we used a polychromatic 
phase retrieval method (see the Methods section) to digitally separate 
the superimposed spectral characteristics of the object holograms. As 
shown in the Results section, this method partially relaxes the tem- 
poral coherence limitation without hardware restrictions and enables 
the use of an illumination source with a wider spectral bandwidth. 
These permit us to work with shorter image acquisition times (i.e., 
higher frame rates) as well as longer object-to-sensor distances. 

Several other computational methods toward coherent imaging 
with broadband X-ray and UV sources have been developed. For 
instance, Parsons et al. treated broadband intensity measurements 



as coherent but noisy images 44 . Furthermore, similar polychromatic 
phase retrieval methods were also reported earlier for lensfree 
imaging using broadband X-ray sources which demonstrated 60 fold 
reduction in the integration time 45 and more than an order of mag- 
nitude increase in usable flux 46 . Our imaging method is different 
from these earlier works in several ways. First, the wavelength of 
operation in our platform is in the visible part of the spectrum and 
our on-chip imaging geometry provides unit fringe magnification for 
the inline holograms of the objects, making our approach signifi- 
cantly wider field-of-view, easily reaching e.g., 20-30 mm 2 or 
>1 cm 2 using CMOS or CCD image sensors, respectively. 
Furthermore, in our reconstruction approach we do not assume 
spectrally invariant objects, and recover a multispectral image cube 
for each object of interest without making small angle or far- field 
assumptions. Finally, each wave propagation step in our reconstruc- 
tion algorithm is performed by using the angular spectrum method 
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Figure 6 | Block diagram explaining our decoding method for lensfree 
multicolor fluorescence imaging. Lensfree, simultaneous multicolor 
fluorescence imaging on a chip based on sparse signal recovery is 
illustrated. 

calculated separately for each instead of rescaling and interpol- 
ating the holograms recorded at a central wavelength as implemen- 
ted in Ref. 45 and 46. Also, our approach does not rely on multiple 
intensity measurements at different wavelengths for phase retrieval 
as reported by Witte et al. toward extreme UV lensfree imaging 47 . 

Reconstructing a spectral image cube from a single intensity mea- 
surement is in general a challenging problem, as a result of which 
PolyPR requires significantly more iterations to converge compared 
to regular PR techniques (e.g., —100 vs. 15). We also investigated 
how thresholding and widening/narrowing the spectral bands affect 
the recovery performance of PolyPR and concluded that using nar- 
row bands selected after a threshold is applied to the intensity of the 
spectrum provides satisfactory imaging performance while also sig- 
nificantly reducing the computational load as illustrated in Supp. Fig. 
1. 

Although spectral demultiplexing and phase retrieval for our on- 
chip holograms are addressing different issues, PolyPR is a combined 
method performing both of these tasks within a single recovery pro- 
cess, run for moderately wide illumination bandwidths. However, 
when the source bandwidth is considerably wider, as in the case of 
sunlight, it is preferred to perform these processes separately. For this 
end, instead of the PolyPR approach for wide-band LED illumina- 
tion, we used a two-step recovery method for lensfree holograms 
recorded with the sunlight illumination. Before any diffraction com- 
pensation (i.e., the holographic back-projection or phase retrieval), 
we performed spectral demultiplexing, where we estimated the 
weights of monochromatic holograms from super- resolved, poly- 
chromatic holograms at three color channels (i.e., red, green and 
blue). Gorocs et al. performed a matrix based cross-talk removal 
for color holograms acquired with red, green and blue lasers 48 ; how- 
ever, our problem is significantly more difficult since we are dealing 
with sunlight that is much more complex and broad in its spectral 
content. We approached this underdetermined problem with a com- 
pressive decoding method by using the knowledge of the combined 
spectrum of the sunlight and sensor sensitivity (i.e., overall pixel 
responsivity) since multispectral spatial data of physical objects are 
naturally redundant 16 ' 49 . 

Another dimension of this work includes simultaneous multicolor 
fluorescent microscopy on a chip. For this end, we demonstrated a 
multicolor fluorescent imaging platform using a custom-designed 
compressive decoding algorithm to recover both the color and loca- 
tion information of fluorescent emitters from their lensfree images. 
In this recovery process, object sparsity is particularly effective as an 
additional constraint since the object distribution (in both space and 
color domains) is already sparse without the need for any additional 
domain representation or transformation. However this process can 
be easily extended into various imaging applications, since almost all 



of the natural and manmade objects of interest can be sparsely repre- 
sented in an appropriate basis 37 with some prior knowledge about the 
objects. Therefore this recovery approach can be applied to other 
fluorescent objects without a technical limitation. 

To demonstrate its proof of concept, we reported in this manu- 
script the imaging results for two different colors; however in prin- 
ciple this method can be extended to fluorescent emitters having 
arbitrary emission spectra within the visible range by using the red, 
green and blue PSFs. On the other hand, due to free-space diffraction, 
the detection SNR is one of the main limitations for extending the 
multi-color imaging performance of lensfree fluorescent imaging on 
a chip. One can in general mitigate this SNR limitation by increasing 
the excitation power level, which could provide a solution using high 
optical density color filters that are installed between the sample and 
the sensor-chip. Note also that most thin-film interference based 
excitation filters or emission spike filters would not work as desired 
within an on-chip imaging configuration since the fluorescent emis- 
sion is not collimated in our lensfree design which would reduce the 
effective optical density of these filters. Therefore, to achieve better 
SNR and contrast in multi-color lensfree fluorescent on-chip 
imaging systems, absorption filters that are customized by coating 
the sensor-chip cover-glass or active area can be utilized to improve 
the multi-color reconstruction performance of the presented 
approach to resolve more complicated multi-color fluorescent sam- 
ples. Another alternative could be to use multi-color quantum-dots 
for fluorescent labeling of specimen such that a deep UV source can 
be used for simultaneous excitation of all the colors in the specimen 
while the sensor- chip remains blind for the excitation wavelength, 
avoiding the use of color filters, also simplifying the design of the 
multi-color lensfree on-chip imager. 

Methods 

Spectral demultiplexing of lensfree on-chip holograms recorded with broadband 
LED illumination. Our polychromatic phase retrieval algorithm (see Fig. 3, top and 
Supp. Fig. 2) iteratively converges to the complex multispectral field distribution by 
using a pixel super-resolved polychromatic intensity measurement (i.e., the pixel 
super-resolved lensfree hologram) along with the source and detector spectral 
characteristics that are known a priori. It starts with the pixel super-resolved 
hologram of the objects 29 , which can be considered as a weighted superposition of 
various monochromatic in-line holograms as detailed earlier. This initial pixel super- 
resolved hologram is scaled with the wavelength weighting factor (c,-) and the square 
root of each scaled hologram is propagated (separately for each A* of interest) to the 
object plane, to calculate the initial estimate. Then, the iterative process is started by 
enforcing an additional constraint, which can vary depending on the nature of the 
object such as phase only, amplitude only, or spatially confined (defining an object- 
support) 50 . Since our objects in this work are both sparse and spatially confined, we 
preferred to use at this step of the algorithm an object-support constraint. For this, we 
use the initial propagation of the hologram intensity at the central illumination 
wavelength to estimate the rough physical boundaries of the object by thresholding 
the resulting field amplitude at the object plane to create a binary object mask. After 
this object-support constraint is enforced, each spatially filtered field is forward 
propagated toward the detector plane, where intensities of each hologram at different 
wavelengths (ki) are appropriately weighted and summed up, and this summation is 
compared to the measured super- resolved hologram intensity by calculating the 
mean square error between the two. At this stage, the phase of the complex field for 
each A-i is kept and its amplitude is updated for each wavelength with the ratio of the 
square root of the calculated and detected intensities. Next, the resulting complex field 
is back propagated to the object plane, completing one iteration, performed separately 
for each A* (see Fig. 3 and Supp. Fig. 2). These successive iterations continue up to e.g., 
100 cycles or simply until the change of the mean square error between the measured 
and calculated hologram intensities is smaller than a certain threshold value, 
whichever comes first. Propagations in both directions are numerically implemented 
by using the angular spectrum method with wavelength specific transfer functions, 
calculated for each Aj 51,52 . 

Spectral demultiplexing of lensfree on-chip holograms recorded with sunlight 
illumination. In our compact and light-weight partially-coherent holographic 
imaging design that is illustrated in Fig. 4, the objects are placed at —0.5 mm away 
from the active region of a color (RGB) sensor array and are uniformly illuminated by 
the sunlight, which is first collected by a simple light collection unit adapted from 
conventional solar cell technologies. After trying some static collection approaches 
such as wide angle lenses or compound parabolic condensers, we decided to use a 
dynamic, angle scanning collection approach (Fig. 4), where we utilized a flexible light 
pipe (1 mm diameter glass core) with a compound parabolic condenser at the 
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Figure 7 | Lensfree simultaneous multicolor fluorescent imaging on a chip. The performance of the decoding method summarized in Fig. 6 is 
demonstrated with (a and b) 4 (am green fluorescent beads and (c and d) 10 (am red and green fluorescent beads, (al, bl, c5 and d5) Raw, lensfree 
fluorescence signatures of these microbeads are cropped from a large field-of-view of —2.42 cm 2 . These raw images are decoded to retrieve color and 
spatial distribution information. (a3, b3, c6, c7, c8, d6, d7, and d8) Decoding results are verified with (a4 and b4 insets) bright-field and (c3, c4, d3, and d4) 
fluorescent microscope images of the same objects. (a2, b2, cl, and dl) The debayered lensfree images illustrate the limitations of the color sensor for 
lensfree fluorescent imaging without our demultiplexing method. The arrows are highlighting two regions, where our method effectively resolves the 
color as well as spatial overlap of lensfree fluorescent signatures. 
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collection end to maximize the light collection efficiency. There is also a plastic 
diffuser after the light pipe to prevent imaging the sky like a pinhole camera would do. 
After the sunlight passes through a plastic diffuser, it is spatially filtered by a 0.1 mm 
diameter aperture, which is positioned at —7.5 cm away from the objects giving us a 
sufficiently large spatial coherence diameter at the sample plane. This 0. 1 mm pinhole 
is controlled using a simple screw-based x-y translation stage, which permits 
capturing slightly shifted in-line holograms of the object. These lensfree images 
captured under sunlight are digitally combined using a pixel super-resolution 
algorithm 31 to form a finely sampled in-line hologram per color channel (i.e., at red, 
green and blue channels). These pixel-super-resolved holograms at each color 
channel are then processed with the knowledge of the spectral response curves of each 
pixel of the sensor array as well as the spectral content of the sunlight to digitally 
remove the cross-talk between color channels and retrieve the multispectral spatial 
features of the objects. Specifically, we minimized the squared 1 2 norm of the 
mismatch between the measurement and the estimation, regularized with the total 
variation of the multispectral object 38 . Then the demultiplexed holograms are 
individually processed either by using iterative monochromatic phase retrieval or a 
simple back-projection if the twin image noise is weak. This recovery process is 
summarized in Fig. 5 (top). 

Spectral demultiplexing for lensfree multicolor fluorescent imaging on a chip. In 

our on-chip fluorescent imaging setup (Fig. l.b and c), the excitation light (center 
wavelength of 480 nm with 15 nm bandwidth) is delivered to the object chamber 
through the side facet of a prism such that after excitation a significant portion of the 
pump beam can be rejected through total-internal-reflection (TIR) occurring 
between the sample holder and the air interface at the bottom. The object chamber is 
placed on the top of an absorption filter (used to reject the scattered excitation 
photons), which makes the object-to-detector distance around 80-100 um. In this 
on-chip fluorescent imaging design, we preferred to use a color sensor (from Kodak, 
KAF-8300) with an active area (i.e., object field-of-view) of 2.42 cm 2 and a pixel pitch 
of 5.4 um. We first experimentally characterized the red, green and blue point spread 
functions (PSFs) of this lensfree imaging system with small fluorescent micro- 
particles (~2 um in diameter). The measured fluorescent intensities from multiple 
isolated particles are interpolated, aligned, and averaged. These intensity profiles 
depend on the structure and optical design of the sensor pixel, filter response and the 
distance between the object and the detector planes. The PSFs for each color channel 
and the Bayer pattern arrangement of the sensor chip are used to decode the color and 
the spatial distribution of the fluorescent emitters from a single raw image by using a 
similar optimization method reported in Ref. 40. However, in this current 
implementation, the forward model is modified to take into account (i) different PSFs 
for different colors (red, green and blue), (ii) under- sampling due to the large pixel 
size, and (iii) the Bayer pattern type, in other words, the specific spatial distribution of 
the color filters on the sensor-array pixels. Figure 6 summarizes this spectral 
demultiplexing recovery process for our lensfree multicolor fluorescent imaging 
platform. 
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